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FOREWORD 


This research was carried out under the partial fin- 
ancial support of NASA through Research Grant NGL 14-005-140 
entitled "Fluid Dynamic Problems associated with Air-Breathing 
Propulsive Systems." Mr. B. H. Anderson and Dr. L. J. Bober 
of the Wind Tunnel and Flight Test Division, Lewis Research 
Center, NASA, served as technical monitors of this research 
grant. This report deals with the problem of a transonic 
flow past a boattailed afterbody under a small angle of 


attack. 



ABSTRACT 


A transonic flow past a boattailed afterbody under a 
small angle of attack has been examined. It is known that 
the viscous effect offers significant modifications of the 
pressure distribution on the afterbody. Thus, the formu- 
lation for the inviscid flow is based on the consideration 
of a flow past a non-axisymmetric body. The full three- 
dimensional potential equation is solved through numerical 
relaxation, and quasi-axisymmetric boundary layer calcu- 
lations are performed to estimate the displacement effect. 

It has been observed again that the viscous effects are not 
negligible. The trend of the final results agreed well with 
the experimental data. It is believed that the general 
three-dimensional turbulent boundary layer calculations must 
be employed before excellent agreement with the experimental 
data can be achieved. 


INTRODUCTION 
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In the attempt to evaluate the performance of an air- 
breathing propulsive system of a supersonic aircraft operat- 
ing at subsonic cruise conditions, an effort has been made 
to calculate the flow field associated with a transonic flow 
past a boattailed afterbody [1,2]* at zero angle of attack. 
Since this afterbody is always immersed in a relatively 
thick boundary layer, built-up along the body from the 
leading edge, it has been found that the pressure distri- 
bution on the boattailed afterbody is always strongly in- 
fluenced by the existence of the viscous layer. It has been 
also learned that the small disturbance potential equation 
is not adequate under this situation to calculate the corre- 
sponding inviscid flow within the transonic flow regime. In- 
stead, the full potential equation must be used and the 
influence of the viscous layer must be accounted for by 
correcting the inviscid flow configuration through the 
displacement considerations. The pressure: distribution on 
the body produced from these calculations showed excellent 
agreement with the experimental data when the flow has not 
been separated and good qualitative indication of the flow 
when it is separated from the rear portion of the afterbody. 
This series cf investigations also fully illustrated the im- 
portance of the interaction between the viscous layer and 
the inviscid flow. 

The present effort, a sequel to this series of investi- 
gations, extends these considerations into the situation 
with small angle of attack. The same axisymmetric model 


♦Numbers in brackets refer to entries in REFERENCES. 
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with boattailed afterbody [1,2,3] is adopted for this study 
whose geometric configuration is shown in Fig. 1. It should 
be noted that this extension leads immediately into con- 
siderable trouble and difficulty. For an axisymmetric body, 
cases with angle of attack are of the three-dimensional 
nature since the flow field also varies in the meridional 
direction. The computational effort (computer time) multi- 
plies with the number of meridional planes under consider- 
ation. In addition, the boundary layer growth along the 
body will also be of a three-dimensional nature, and the 
corrections to the inviscid body configuration would lead 
itself into a non-axisymmetric one even if the solid body is 
of the axisymmetric shape initially. 

For an axisymmetric body, the coordinates to describe 
such a flow problem are naturally of the cylindrical system. 
Due to the flow with angle of attack, the radial velocity 
component along the axis will no longer vanish and some of 
the terms in the general three-dimensional potential equation 
in this coordinate system will approach infinity.* It is 
thus obvious that this equation cannot be employed for the 
flow along the axis and additional manipulations are needed 
to overcome this problem. 

Since the compressible turbulent three-dimensional 
boundary layer calculations have not been developed, it is 
assumed that under the condition of small angle of attack, 

*It is, nevertheless, interesting to observe that the uni- 
form approaching flow condition with angle of attack satis- 
fies this mathematical expression describing the potential 
equation in the cylindrical system of coordinates. 
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the viscous portion of the cross flow may not be very im- 
portant, and the boundary layer growth can be approximated 
by a quasi-axisymmetric viscous flow calculations. Thus, 
the main purpose of the present effort is to study the 
feasibility of carrying out the three-dimensional transonic 
flow calculations, and developing the computer program for 
its inviscid flow. It is felt that it is simple to replace 
in the future the quasi-axisymmtric viscous flow treatment 
by a full three-dimensional compressible turbulent boundary 
layer calculation when it becomes available. Indeed, it has 
been found from the results of the present investigation 
that the equivalent inviscid configuration of the afore- 
mentioned boattailed afterbody is sensitive to the influence 
of the viscous flow, and accurate results can only be 
produced from the three-dimensional turbulent boundary layer 
considerations. The detailed treatment of this problem in 
its present form is described in the following sections. 


I 
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BASIC GOVERNING EQUATIONS 

Inviscid Flow 

Referring to Fig. 2 where an axisymmetric body at an 
angle of attack is depicted, the compressible three-dimensional 
potential equation in the cylindrical system of coordinates 
(z,r,w) governing the flow is given by 

II - (* 2 /c 2 )) 4 ZZ + U - (» 2 /c 2 )] J rr + [1 - (t 2 /r 2 c 2 )] 

• ( W r2) - I< 2 W /(r2 c2)) 

- 4 z )/(r 2 c 2 )] - £(2* r t z )/c 2 ] 

• * rz + <» r /r) 11 + (t 2 /r 2 c 2 )J = 0 (1) 

where the subscripts denote partial differentiations and c 
denotes the local speed of sound which is related to the 
stagnation velocity of sound c q and the free stream flow 
condition through 

c 2 = cl - [(y - l)/2][$ 2 + <t» 2 + $ 2 /r 2 ] (2) 

O a Cl) 

2 2 2 

with c q = c^ + [ (y - l)/2] v^. Upon introducing 

4> = v^ (z cos a + r cos to sin a + <$>) , the disturbed potential 

<j> is found to satisfy 

[1 - (U 2 /C 2 )] 4> zz + [1 - (V 2 /C 2 )] <f rr + [1 - (W 2 /C 2 )] 

* I (4» - r cos to sin a)/r 2 ] - (2VW/C 2 ) 

* l (<J> - sin to sin a)/r] - (2UW/C 2 ) (<f> /r) 

XTtJu CO c* 

- (2UV/C 2 ) <J> zr + (V/r) [1 + (W 2 /C 2 )] = 0 (3) 


where 
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u a a cos a + <|> z ; 

V = 4> /V = cos w sin ct + <i> ; and (4) 

W “ $ /r v = 6 /r - sin w sin a 

(A) 00 (A) 

and 

C 2 = (c 2 /V 2 ) = 1/M 2 + I(y - l)/2) (1 - U 2 - V 2 - W 2 ) 

In anticipating that the viscous effect will offer im- 
portant corrections to the configuration of the body, and 
this correction varies in the meredional direction, ono thus 
expects that the effective body is non-axisymmetric, namely, 
r b is a function of both z and w. One may introduce trans- 
formations according to 

£ = z ; x\ - r - r b ( z , w ) ; and Q = <d (5) 

and an additional transformation of the form 

S = [Bn/(1 + Bn)] (6) 

so that the condition of n ® is equivalent to £ = 1 in the 

finite computational space. B (B = 8) is a stretch parameter. 
Equation (3) appears now in the new system of coordinates 
(£,S,G) as 

DD<f>££ + EE* B 2 (1 - ;) 4 + FF<j) fifi 

+ B(i - + B (1 - C) 2 *EF<J>^ + FD<^ 

+ B (1 - C) 2 {E - 2B (1 - 0*EE] <j> 5 + G = 0 (7) 

where 

DD = [1 - (U 2 /C 2 )] (8a) 

EE = [I - (U 2 /C 2 )] r 2 z + [1 - (V 2 /C 2 )] 

+ [l - (W 2 /C 2 )] (rj^/r 2 ) + [(2VW r bu /(r C 2 )] 

- K2WU r bz r b J/(r C 2 )] 

+ (2UV/C 2 ) r bz 


(8b) 
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FF « [1 - (W 2 /C 2 )]/r 2 (8c) 

DE - [2WU/(r C 2 )] r bw - 2[1 - (U 2 /C 2 )] r b2 

- (2UV/C 2 ) (8d) 

EP * [2UW/(r C 2 )) r b2 - [2VW/(r c 2 )] 

- (2(1 - W 2 /C 2 ) r bw ]/r 2 (8e) 

FD = - [2WU/(r C 2 )] (8f) 

E » [2UW/ (r C 2 )J r ba)2 - [1 - (U 2 /C 2 )] r b22 

- [1 - (W 2 /C 2 ) ] (r bUM /r 2 ) (8g) 

G = [1 + (W 2 /C 2 )] (V/r) 1 [2VW/(r C 2 )] sin u sin a 

- [1 - (w 2 /C 2 )] [(cos a) sin a)/r] (8h) 


with r = r b + ?/[3(l - £)] and the subscripts z,w for r b 
indicate respective partial differentiations of r^. In the 
transformed coordinate system, U, V, and W are given by 

U = cos a + <J>£ - r b2 B(1 - c) 2 ^ (9a) 

V = cos ft sin a + B (1 - c) 2 (9b) 

w = d/r) t<f n - r bu B (1 - C) 2 * ; ] 

- sin ft sin a (9c) 

The boundary condition for the flow can be simply stated as 

V * n = 0 on the body surface where n is the unit outward 

normal. Referring to Fig. 3, this condition of tanger.'-.v on 

the surface can be further developed as follows: 

The unit normal n can be expressed as 

n = -cos 8’ sin 0 e + cos 8' cos 0 e + sin 8* e 

z r u) 
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where tan 6 = (dr^/dz). The boundary condition thus can be 
expressed as 

V ■ U tan 0 - W tan fi'/cos 6 at ; * 0 

since 

(tan P'/cos e) ■ tan B - - (l/r b ) (8r b /3w) » -(1/^) r bw 

and upon substituting Eq. (9) into the above expression, the 

boundary condition becomes 

4^(0) = I(4>^(0) + cos a) r faz + (<f> fi (0)/r b “ sin sin 
• (l/r b ) r ba) - cos n sin a]/(B[l + r 2 z + (1/r 2 ) r 2 u ] } (10) 


where (0) signifies that the quantities are to be evaluated 
at ; = 0. 


The boundary condition far awr.y from the body now 
becomes 


<J> = 0 at ; = 1 or £ = +<° (11) 

However, in view of the experience gained from previous 
study [1] on the flow without angle of attack, it is appro- 
priate to assume that $ = 0 at £ = -2 and £ = + 3 where the 
length, L, of the body (Fig. 1) has been chosen as unity. 

Equation (7) may now be solved by finite difference 
calculdtic-.jj through relaxational method. However, as i en- 
tioned in INTRODUCTION, Eq. (7) cannot be applied for the 
points along the axis anead of the body. The governing 
equation in the Cartesian system of coordinates must be 
introduced for all the points on the axis. This equation is 
given by 


[1 - ($*/c 2 )] * z2 + [1 - ($ 2 /c 2 )] 4> xx 


. * - [(2* $ )/c ] $ - 1(2$, 


yy 


XZ 


- 1(2* V /C 1 *yz 


+ H ■ (^y/C 2 )] 

t )/c 2 ] $ 
y " J xy 


0 


( 12 ) 


i 


Since this equation is only applied to points on the Z-axis 
ahead of the body and the fact that the flow is symmetric 
with respect to the x-z plane leads to ♦ ® 0 for all points 

on the axis, Eq. (12) is immediately simplified to 

U - <^/= 2 » ♦„ + 11 - <» 2 /c 2 )J » xx ♦ * yy 

- I(2* x y/c 2 ] 4> xz * 0 (13) 

The disturbed potential function, <p, is introduced according 
to 4> = V w (z cos a + x sin a + <f>) which is precisely the same 
disturbed potential function <p introduced previously. 
Equation (13) becomes 

{1 - [(cos a + $ Z ) 2 ]/C 2 } <P zz + {1 - [(sin a + 4> X ) 2 )/C 2 } 

* *xx + ^yy " {[2 (cos a + Q g ) (sin a + 4> x ) ] /C 2 > x2 = 0 

(14) 

where 

C 2 = (1/M 2 ) + [ (y - l)/2] [1 - (cos a + y 2 
- (sina+ <(> x ) 2 ] 

When Eq. (14) is written in finite difference form with the 
central difference scheme, the grid points involved for a 
typical point on the z axis are shown in Fig. 4. It is 
recognized that all grid points are already the points 
needed in the finite difference formulations of Eq. (7) ir. 
the transformed system of coordinates (£*£«(?)• No ad- 
ditional new grid points are introduced for the calculation 
of Eq. (14) . 
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Finite Difference Equations 

The computational space of -2 ^ + 3, 0 <_ c £ 1, and 

0 < ft < 7T was discretized into a grid system where uniform 
values of As * 0.04 and AQ = 30° = 0.523 rad were employed. 
Thus, calculations on seven meridional planes have been 
carried out. Since the present investigation has the main 
interest in the flow properties within the boattailed after- 
body region, variable values of AC were used. The grid 
system is shown in Fig. 5. 

Many of the finite difference forms of the terms in Eq. 
(7) are different for locally subsonic or supersonic flow. 
They are given by 

“ {ACj/U^AC-l + ac 2 )]} *J +1#j 

+ I(AC 2 - AC 1 )/(AC 1 A£ 2 )] 4^ 

- {AC 2 /[AC 1 (AC 1 + AC 2 )J} ^_ 1#j 
= AB1 ^ +lfj + AB2 + AB3 ^ . (15a) 

4^ = {2/[AC 2 (AC 1 + AC 2 )]} -U2/(AC 1 AC 2 )] )* 1 l t ■ 

+ {2/(AC 1 (AC 1 + AC 2 )H j 

“ BB1 +i + i.j + BB2 + BB3 ♦ti.j 


(15b) 
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♦«- <i/ 2 «> !“!(♦][«, j- ♦i+i.j* 

+AB2(<{>* + ^ - **“*) + AB3(<j>£J . - <j£ 

l r J 1, J 1 ± f J 1 

♦CD - U/2AG)lABl(** +lij+1 - ♦J +Xfj . 1 ) 

+ AB2 <j + l - *2,j-l> 

+ AB 3 (*i. Xi j+i " 


for locally subsonic flow, and 


♦ 5 = {[AC q + 2AC 1 ]/tA^ 1 (AC 0 + AC X )]} 4>J fj 


- l(A5 1 + A5 0 )/(A5j. A? 0 >] *2-1,3 


+ <IAt 1 )/U5 0 (A5 0 + 45x11) *2_ 2>j 


API *2 (j + *P2 *2-1, j + AP3 *2-2,j 


* ?5 = { 2/ [45x (45 0 + 45x11 } 


- I2/(45 0 45x)1 *2-1,3 


{2/[45 0 (45 0 + 45x11} *2- 2 ,3 


BP1 *2 rj ♦ BP2 *2_x, j + BP3 *2. 2jJ 


j ,)] (15c) 

1 • J 


(15d) 


(16a) 


(16b) 
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♦« - (1/2AC)IAP1(^ - ♦£}) + *P2U5!i tJ - ♦ili.j' 


+ " 3 ‘*i-a.J ' ♦Ea.j” 


(16c) 


♦{0 = U/2ABHAPl{*5 (j+1 - *5^.!) 

+ AP2 <*i-l,j + l - 
+ AP3($j_ 2j j +1 - 

for locally supersonic flow where 

A ^2 = ? i+l " ^i ; 1 ^i-1 ' and 

A ^o = 4 -1 “ ^i-2 
Other finite difference forms are 


(16d) 


(17) 


<f>. = A; 

<? 1 r J / J 

6 = (<f> k+1 - 2(j) k . + <{) k-1 ) / (Ac) 2 

vv i,D v i#D y i/D 

♦a ‘ <i + i - "'i.j-i’'' 2 "’ 


(18a) 

(18b) 

(18c) 



<j >4 , ,)/[ (An) 2 ] 

i# J J* 


(18d) 
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♦OB’ ‘♦i.J+l ‘ 2 *i,J 


,k+l 


,k-l 


k+1 


.k-1 




[ (<}>." - <J>? * ) - U. t . - <h . t .)]/ 

Y i/3+l y i,D+ 1 ' v i,;j-l y x,3~l 


[4AC Aft] 


(18e) 


Equation (7) may now be written in finite difference form 
for the point (£.,io. , c v ) and the resulting expression, when 
arranged in tridiagonal form for all points along the c co- 
ordinate with fixed values of is 

[B 2 (1 -;) 4 *EE + (Ac/2) B (1 - ;) 2 *DE • AB2] 


* + IAS 2 • DD • BB2 - 2B 2 (1 - ;) 4 • EE 

- 2 (Ac;/ AR) 2 • FF] & . + [B 2 (l - ^ ) 4 • EE 

J 

- (AC/2) B (1 - C) 2 • DE • AB2] = 

J 

- (AC/2) B (1 - C) 2 * DE • [ABl j - 

+ AB3(<J>*+^ . - - (AC) 2 

• (DD • (BBl + BB3 (j)*.^.) + (FF/ (Aft) 2 ] 


‘♦i.j+i + + “) 


(19) 
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for subsonic flow, and 

[B 2 (l - 0 4 • EE + (Ac/2) B (1 - 0 2 * DE • API] 

• <jf + J + I(Ac) 2 • DD • BPl - 2B 2 (1 - ;) 4 • EE 

X 9 J 

- 2 (A;/ Aft) 2 • FF] <jK . + IB 2 (1 - c) 4 • EE 

i# 3 

- ( AC/2) B(1 - O 2 • DE • API] = 

■W J 

- (A;/2) B(1 - C) 2 • DE . [AP2(4,^ . - 

1” J- / J 

.k-1 v . k+1 ,k-l 

+ AP3( +i-2,j ' 

- (At) 2 IDD • (BP2 ^ 1 1 + BP3 *£ , ..) 

X“ i f J X-4, t J 

+ Iff/ ( A ft) 2 ] + (20) 

for supersonic flow where 

H = B(1 - f) 2 • EF • + a? + FD • * a5 
+ B (1 - ;) 2 IE - 2B (1 - C) • EE] $ + G (21) 

To account for the boundary condition on the body 

surface (c = 0,k = 1), 4^(0) was evaluated for the points on 

the surface according to Eq. (10) and a row of grid points 

was introduced at c = -Ac in every meridional plane whose 

values of <j> (identified as ^ . ) were given by 

1 » 3 


( 22 ) 
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*i,j " *i7j " 2Lz > c = 0) 

Similarly, due to symmetry, <j>^ and r^ vanish for all 
grid points on the meridional planes of u = 0 (j=l) and 
u = n(j = 7). Again, a row of grid points was introduced 
for each 5 value at u = -Au or w = it + Aw whose values of 

k 

<|> (identified as <J>. . ) were given by 

1 , D 

*i,b ■ *i,j =2 £or “ - 0 < 23a > 

K - * for 0 ) = TT (23b) 

When Eq. (19) or (20) is written for all points along the 
5 coordinate (for fixed values), the coefficient 

matrix for <J> is tridiagonal, and values of <j> can be solved 
by standard and efficient methods. Line relaxational cal- 
culations can be carried out along each plane of fixed £ 
value in the direction of increasing u and then sweeping 
from upstream toward downstream positions, and the information 
of the fictitious row of grid points created to account for 
the boundary condition was updated as soon as the calculations 
for that column of grid points were completed. However, 
before carrying out these calculations, finite difference 
calculations must be performed for the points along the z- 
axis ahead of the body. 

Referring to Fig. 4, the finite difference forms of the 
terms in Eq. (14) are given by* 

*The finite difference forms are for subsonic flow ahead 

of the body. For supersonic flows, expressions for $ , 

<J> , and d> must be modified. 

T zz' Y zx 
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♦ = ABl ^? +1 + AB2 <()P + AB3 4>P_ 1 (24a) 

4> ■ BB1 (J)?. , + BB2 4>J + BB3 4>? - (24b) 

ZZ T l+1 1 T 1 -I 

♦zx = U/(2A ni )] [AB1(^ +1 - +® +1 ) 

+ AB2(<j>" - <f>®) + AB3 (<j ) ^_ 1 - 1 (24c) 

4>x = (♦? " ♦J')/2An 1 (24d) 

4» xx = (4>" - 24 .P + *J)/[(An 1 ) 2 ] (24e) 

<J>yy - 12(4>® - ^))/[(A ni ) 2 ] ( 2 4 f ) 


where Ar^ corresponds to A; between the first and the 
second grid points along the 5 coordinate and superscripts 
p,n,s, and e refer to the points on the axis, north, south, 
and east of the axis point, respectively. Upon writing Eq. 
(14) in finite difference form, the value of 4>? can be 
readily obtained from 
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<J>? = KAn-j^ 2 {(2 (cos a + 4> z ) (sin a + <j> x )/C 2 ] <f> zx 

- 11 - (cos a + 4> z ) 2 /C 2 } [BB1 <J>P +1 + BB3 

- U - (sin a + <P X ) 2 /C 2 : (4-J + <pf) - 2<f>J]/ 
[(An 1 ) 2 BB2 [1 - (cos a + <j> z ) 2 /C 2 ] 

- 2(1 - (sin a + <p x ) 2 / C 2 ] - 2] (25) 


Viscous Flow 

The quasi-axisymmetric boundary layer calculations are 
based on the integral method of Sasraan and Cresci [4] and 
follow in the same formulations as described in (1] . The 
pair of differential equations for the momentum thickness 0 c 
and the form factor H c are integrated numerically with the 
given inviscid surface flow condition along the axial direc- 
tion as the guiding free stream. These equations are pre- 
sented again as follows: 

(df/dx) = 1.268 [ (-f/M ) (dM /dx) (1 + g H. ) 

c C W X 

- (f/r b ) (dr b /dx) + A] (26) 


and 
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dH./dx » -(1/2M ) (dM /dx) [H. (H. + l) 2 (H. -1)] 

1 6 6 JL A 

• U + (g w “ 1)[(H 2 + 4H i - l)]/[(H i + 1) (H^ + 3)]} 
+ [(H 2 - l)/f] A{H i - [0.011(H i + 1)0^ - l) 2 / 
H 2 ] • (2/c f ) (T 0#e /T)} (27) 


where 


f = r (M e c Q 6 i )/v o ] 1,268 ; (28a) 

e i = 0(T e //T o e ) (Y +1 )/l 2 (y-D] (28b) 




E T /T 
vr o,e 


(28c) 


A - 0.123 e' 1-561 H 1 


• t< M e c o )/u o J 1T e /fl IT 8' /T o,e l3 ‘^o 1 °' 263 (28d) 
T/T 0>e = 0.5 (T„/T 0 _ e ) + 0.22 Pr 1/3 

+ (0.5 - 0.22 Pr 1/3 ) (T /T ) (28e) 

6 Uf c 

H c - H. g w Il + ( <Y - l)/2] M*J 

+ [ (Y - 1 >/23 ( 28f ) 

I 
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(c f /2) = 0.123 e“ 1,561 H i [M e c Q O./vj" 0 * 268 (T g /T) (p/u o 
and x is the length along the surface of the body. 




J 0 .268 

(28g) 


METHOD OF CALCULATIONS 
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The method of calculations follows in the same manner 
as that described previously for cases with zero angle of 
attack [1] . The successive refinement of grid spacing in 
the ^-direction is again employed and its number of grid 
points changes from 26 to 51 and finally to 101. Initially, 
calculations are carried out for all grid points of the same 
K location (they are all on the plane perpendicular to the 
z axis) by performing line relaxations from u = 0 to u = u , 
followed by all grid points in the next downstream constant 
£ -plane, and so on. The information at the grid point 
created to account for the boundary condition on the surface 
is updated, according to Eq. (22) immediately after the 
calculations for that column of grid points are completed. 

As the end of a complete cycle of relaxation for all £ 
stations, values of for all points on the axis ahead of 
the body are updated according to Eq. (25) . From the ex- 
perience gained in the study of cases with zero angle of 
attack, frequent quasi-axisymmetric boundary layer calcu- 
lations are not required. These calculations were carried 
out at the end of second, fifth, ninth, fourteenth, twenti- 
eth,..., cycle of inviscid flows calculations and the 
corrections to the inviscid body configuration through the 
displacement effect of the viscous layer were applied ac- 
cordingly (see Ref. [1] for details of this correction). After 
the successive grid refinements, the final flow pattern is 

eventually established when the successive change of 4> is 

-5 

less than a small number (c = 3 x 10 )* for all grid points 

throughout the field. 

*To save computational time, this number is ten times larger 

than the valnp emnlnved for raspc in Ref. Ml . 


RESULTS OF CALCULATIONS 

After completion of developing the computer program 

for flows v/ angle of attack, one possible way of checking 

the program s to reproduce the results when the angle of 

attack is rot to zero. Indeed, the results obtained for 

different meridional planes at = 0.8 are agreeable to 

each other and are shown in Fig. 6 where the previous re- 

— 6 

suits with e = 3 x 10 (Ref. {1] ) and the corresponding 
experimental data are also presented for comparison. The 
slight difference between the two calculations appears to 
come from the difference in c for the criterion of con- 
vergence. 

Calculations for the case of flow with angle of attack 
with a = 4° at M =0.8 have also been carried out. The 

CD 

inviscid and viscous flow results are presented in Fig. 7. 
Inviscid flow results were obtained by bypassing the boundary 
layer calculations. It is obvious that the viscous flow 
modifies the pressure distribution on the boattailed after- 
body so much that it is impossible to obtain reasonable 
results without consideration of its effects. The experi- 
mental results obtained by Shrewsbury (5] corresponding to 
this flow case are presented in Fig. 8 along with the vis- 
cous flow results for comparison purposes. Although the 
viscous layer is calculated from the quasi-axisymmetric 
approach, the final predicted level of the pressure coef- 
ficient, as well as its trend of variation is reasonable in 
comparison with experimental data. It is believed, however, 
that the quasi-axisymmetric boundary layer calculations must 
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be replaced by the three-dimensional boundary layer calcu- 
lations before excellent agreement can be achieved. The 
corresponding displacement thickness resulting from the 
quasi-axi symmetric boundary layer calculations are shown in 
Fig. 9. The computational time required for this case is 
225 seconds on the CDC CYBER 175 computer system. Since the 
computational cost is expensive while tne method of viscous 
flow analysis is not accurate, no extensive calculations 
have been performed. 
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CONCLUSION 

A computer program dealing with an axisymmetric body 
at angle of attack within the transonic flow regime has been 
developed. The viscous layer growth on the body introduces 
important modifications to the inviscid body geometry within 
this flow regime. A quasi-axisymmetric boundary layer esti- 
mation serves to show the correct trends of the pressure 
distribution but would not produce accurate results when com- 
pared with the experimental data. A fully three-dimensional 
compressible turbulent boundary layer treatment must be imple- 
mented to improve this situation. Nevertheless# the investi- 
gation does establish the fact that inviscid flow field of an 
axisymmetric body immersed in a transonic flow at angle of 
attack can be successfully calculated. 
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NOMENCLATURE 


A,f 


Afil,AB2, 

AB3 


API, AP2, 
AP3 


B 

BB1 , BB2 , 
BB3 


BP1,BP2, 
BP 3 


C 

C. 


functions for boundary layer calculations, 
defined in Eqs. (26) and (28) 


coefficient functions at point where flow is 
subsonic, defined in Eq. (15a) 


coefficient functions at point where flow is 
supersonic, defined in Eq. (16a) 

radial stretch parameter (B = 8) 


coefficient functions at point where flow is 
subsonic, defined in Eq. (15b) 


coefficient functions at point where flow is 
supersonic, defined in Eq. (16b) 


DO , EE , DE , 
FF,EF,FD, 
E,G,H 


max 

D s 

H c 

H i 

i, j»k 
L 
M 
Pr 

P 

R 


dimensionless speed of sound 
pressure coefficient (p - p^J/d/^) p TO V 
local speed of sound 
skin friction coefficient 


coefficients in inviscid analysis, defined 
Eqs. (8) and (21) 

maximum model diameter 

sting diameter 

shape factor of compressible boundary layer 

shape factor of incompressible boundary layer 

indices for finite different formulations 

body length to end of boattail 

Mach number 

Prandtl number 

pressure 

boattail juncture radius of curvature 
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r 

r b (z,u) 

T 

T 


U,V,W 


u,v,w 


X 

Y 

6 * 

AC 

AC 


n 

0 

U 

y 

v 

C 

* 


(0 

ft 

a 

0f$' 


radial distance from the axis 
local body radius 
temperature 

reference temperature for the boundary- layer 
calculation 

dimensionless velocity components in axial 
direction, radial and meridional directions, 
respectively 

velocity components in axial, ’-adial, and meri- 
donial directions 

free stream velocity 

distance along body surface, or Cartesian coordinate 
Cartesian coordinates 
ratio of specific heats 

displacement thickness of the boundary layer 

uniform spacing of the finite grid system in 
the c-direction 

nonuniform spacing of the finite grid system 
in the ^-direction 

transformed coordinate, r - r^ 

momentum thickness of the compressible and the 
cooresponding incompressible boundary layers 

viscosity 

viscosity evaluated at T temperature 
kinematic viscosity 
transformed axial coordinate 
potential function 

normalized disturbance potential function 
meridional coordinate 
transformed meridional coordinate 
angle of attack 
coordinate angles, see Fig. 3 


Subscripts 


e 


r,z,w 


w 

C 

o 


OP 


edge of boundary layer 

partial derivatives with respect to r,w,z, 
respectively 

wall 

partial derivative with respect to £ 
stagnation state 
approaching free scream 


x,y,z 

partial derivatives with respect to 
x,y,z, c, respectively 


Superscripts 

p,n,s,e grid points on the axis ahead of the body 
north, south and east of this point 


